Preprint typeset in .IHEP stylo - HYPER VERSION 



hep-th/0606086 
HUTP-06/A0022 
MIT-CTP/3754 



Ricci flow and black holes 



Matthew Headrick 

Center for Theoretical Physics, Massachusetts Institute of Technology 
77 Massachusetts Ave., Cambridge MA 02139, USA 
E-mail: jieadrickOmit . edn 



Toby Wiseman 

Jefferson Physical Laboratory, Harvard University 
Cambridge MA 02138, USA 



E-mail: twisemanOf as .harvard, edu 



Abstract: 

Gradient flow in a potential energy (or Euclidean action) landscape provides a natural set of 
paths connecting different saddle points. We apply this method to General Relativity, where 
gradient flow is Ricci flow, and focus on the example of 4-dimensional Euclidean gravity 
with boundary S"^ x S*^, representing the canonical ensemble for gravity in a box. At high 
temperature the action has three saddle points: hot flat space and a large and small black 
hole. Adding a time direction, these also give static 5-dimensional Kaluza-Klein solutions, 
whose potential energy equals the 4-dimensional action. The small black hole has a Gross- 
Perry- Yaffe- type negative mode, and is therefore unstable under Ricci flow. We numerically 
simulate the two flows seeded by this mode, finding that they lead to the large black hole and 
to hot flat space respectively, in the latter case via a topology-changing singularity. In the 
context of string theory these flows are world-sheet renormalization group trajectories. We 
also use them to construct a novel free energy diagram for the canonical ensemble. 
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1. Introduction and summary 

Gradient flow is a powerful theoretical tool for exploring a system's potential energy landscape, 
which has been applied in contexts ranging from biophysics to string theory. Similarly, it is 
sometimes useful to study gradient flow of a system's action on its space of histories, for 
example in the context of the path integral formulation of quantum mechanics. (Like the 
path integral, the gradient flow will strictly speaking be well defined only in the Euclidean 
framework.) 

The purpose of this paper is to investigate gradient flow, both of the potential energy and 
of the Euclidean action, in the context of General Relativity. Gradient flow of the Einstein- 
Hilbert action on the space of metrics on a given manifold is governed by the following 
equation: 

%^ = (1.1) 



- 1 - 



where A is the flow time. In the mathematical literature, the flow defined by equation (^]^) is 
known as Ricci flow, and in recent years has been the subject of intense study. It also arises 
in string theory as the one-loop approximation to the renormalization group (RG) flow of 
sigma models. Section 2 of this paper is devoted to a review of the basic properties of Ricci 
flow. 

As our laboratory for the application of Ricci flow to gravity, we take one of the simplest 
systems that admits multiple saddle points: 4-dimensional pure gravity (without cosmological 
constant) in a spherical box of radius i?, in the canonical ensemble at temperature (3~^. The 
path integral for this system involves all Euclidean metrics, on all topologies, with boundary 
5^ X Sj^. At high temperature there are three saddle points of the action: a large and a small 
black hole, and the product of with a flat 3-ball. The Lichnerowicz operators for the large 
black hole and hot flat space metrics are positive definite, whereas, as we show, the one for 
the small black hole has a single negative eigenvalue, analogous to that of Gross, Perry, and 
Yaffe (GPY) in the asymptotically flat case. This indicates that the small black hole is 
truly a saddle point of the action, whereas the other two solutions are actually local minima 
(not counting Weyl fluctuations of the metric, which always decrease the action; these are 
dealt with in the path integral by rotating their contours of integration). GPY argued that, 
due to this negative mode, the small black hole serves as a bounce (or sphaleron) that allows 
the system to pass from one minimum to the other by thermal fluctuations. Furthermore, 
according to arguments given by Whiting and York g. Whiting g, Prestidge ||, and Reall 

this negative mode is related to the local thermodynamic instability of the corresponding 
Lorentzian black hole. 

We can also add a time direction to the system, and consider 5-dimensional Lorentzian 
gravity with boundary R x 5^ x 5^. Ignoring the Hamiltonian constraint (which in any case 
is automatically satisfied at the saddle points), the space of time-symmetric initial data is the 
space of Euclidean 4-metrics, and the potential energy is the 4-dimensional Euclidean action. 
The three saddle points give static solutions, which are reinterpreted in this context as large 
and small bubbles of nothing and the Kaluza-Klein vacuum, respectively. Their dynamical 
stability is controlled by the Lichnerowicz operators for the respective Euclidean 4-metrics, 
so the small bubble has a tachyon while the other two solutions are stable. Further details 
about both the Euclidean and the Lorentzian versions of this system are given in Section 3. 

In both of the above contexts, Ricci flow describes gradient flow: of the potential energy 
in the 5-dimensional Lorentzian context, and of the action in the 4-dimensional Euclidean one. 
The small black hole's negative mode indicates that it is unstable under Ricci flow. This raises 
the question of where the two flows that start at this saddle point, moving opposite ways along 
the unstable direction, ultimately end up. One guess would be that each flow asymptotes to 
one of the other two saddle points. Using numerical simulations we were able to confirm this 
guess; this is the main result of this paper. Ricci flow thus defines a unique path in the space 
of metrics connecting the small black hole to each of the other two saddle points. The one 
connecting the small black hole to hot flat space is particularly interesting, because at a flnite 
flow time the metric passes through a localized topology-changing singularity (the metrics 
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are otherwise smooth throughout the flows). The flows are described in detail in Section 4. 

These flows have various applications. In the context of string theory, they represent 
world-sheet RG trajectories, with the small black hole as the ultraviolet fixed point and hot 
flat space and the large black hole as the infrared ones. It is interesting that they qualitatively 
mirror the time evolution of the Lorentzian small black hole perturbed by its thermodynamic 
instability, as it either evaporates or grows to become the large black hole. Quantitatively, 
however, the evolution is different, since at every point along the evaporation process the 
geometry remains Schwarzschild to a good approximation, whereas as we will see along the 
RG flows it does not. 

RG flow is also often applied to the spatial part of a world-sheet theory as a technique 
for studying closed-string tachyons. Hence we may use the terminology of the 5-dimensional 
Lorentzian context and say that the small bubble of nothing (which has a tachyon, and there- 
fore a relevant operator on the world-sheet) is the UV fixed point, and the large bubble and 
Kaluza-Klein vacuum are the IR fixed points. These flows are in accord with the arguments 
made in the paper j^, that the ADM energy of the target space should be lower in the IR than 
in the UV of a world-sheet RG flow. The full dynamical evolution of the unstable small bubble 
has been performed by Lehner and Sarbach j^, |8|, and we may contrast this behavior with the 
RG flows. In the case where the small bubble initially expands, the dynamical simulations 
were consistent with the large stable bubble being the end state, although the simulations 
were not performed in a box and so saw only continued expansion. The more interesting 
case is initial collapse of the small bubble, which in RG flow goes to Kaluza-Klein flat space, 
but dynamically forms a Kaluza-Klein black string. Whilst this will indeed evaporate to flat 
space eventually, RG flow and classical dynamics clearly do not agree for this system. We 
attribute the difference to the fact that this flow involves a topology change, and whilst this 
may occur smoothly in RG flow, in the dynamical case cosmic censorship is expected to form 
a horizon to shield it. 

Another application of our flows is in the Euclidean (thermal) context. The fact that 
Ricci flow is gradient flow means that, at any point along the flow, the action is stationary 
in the directions orthogonal to it. Hence in principle these directions can be integrated out 
to define an off-shell free energy; in a saddle point approximation this free energy is simply 
equal to the value of the action at that point on the curve. Thus, Ricci flow can be used to 
construct a new type of free energy diagram for this thermal system. This diagram, shown 
in figure 10, is discussed further in Section 5. 

Although we studied the specific case of 4 dimensions, it seems likely that the flows 
are qualitatively similar in higher dimensions. Also, while we employed a simple box as 
our infrared regulator, our results would likely have been qualitatively unchanged had we 
instead included a negative cosmological constant and imposed asymptotically anti-de Sitter 
boundary conditions. It would be interesting to see whether they have physical implications in 
the context of the AdS-CFT correspondence, and we conclude the paper with a brief comment 
on this. 

Ricci flow has previously been studied numerically in a variety of contexts. For example. 
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a 2-dimensional flow asymptoting to the dilaton black hole was studied by Hori and Kapustin 
The formation of singularities in 3 dimensions was studied by Garfinkle and Isenberg 
1^, 0]. And Ricci flow was investigated as a numerical method for solving the Einstein 
equation on the 4-dimensional manifold K3 by the present authors [12|. However, as far as 
we know the flows we study in this paper are the first explicitly known that connect different 
saddle points of the action. 



2. Ricci flow 

Ricci flow is an analog of the heat equation for geometry, a diffusive process acting on the 
metric of a Riemannian manifold. In this section we define it, explain its basic mathematical 
properties, and describe how it appears in several different contexts in physics. Indeed, histor- 



ically it was in physics that Ricci flow first made its appearance, in the work of Friedan |13| on 
the renormalization group for two-dimensional sigma models. It was subsequently introduced 
into mathematics by Hamilton as a tool for studying the topology of manifolds. Among 
other things, he developed a program, recently brought to fruition by Perelman [0, for 
proving Thurston's geometrization conjecture concerning the classification of 3-manifolds (for 
reviews see |l^, |l^). By now there is a very extensive literature (including a textbook |p!9| ) 
on Ricci flow. 

The material in this section is review, except the discussion of boundary conditions in 
subsection 2.2, which is new as far as we know. 

2.1 Definition and basic properties 

Let M be a D-dimensional compact manifold, and Ti. the space of Euclidean (i.e. positive 
definite) metrics g^u{x) on it {not the space of metrics modulo diffeomorphisms) . Ricci flow 
is the first-order flow on with respect to an auxiliary "time" variable A, defined by 



dX 



-2R^,ix,X). (2.1) 



The flow is obviously invariant under (A-independent) diffeomorphisms of M. It is sometimes 
useful to consider the more general equation 



dX 



-2i?^, + 2V(^e.), (2.2) 



where ^'^{x, A) is some vector field. In fact, it's easy to see that the solutions to equations (|2.lD 
and (|2.2[ ) are the same, up to the diffeomorphism obtained by integrating with respect to 
A. Thus the two equations describe different flows in 7i, but the same flow in the quotient of 
7i by diffeomorphisms of M. 

In what sense is Ricci flow diffusive? Let us examine how a short-wavelength, small- 
amplitude perturbation of the metric evolves. If the wavelength is much shorter than the 
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typical curvature radius of the manifold, then it is sufficient to expand about flat space. With 
9fj,i/ = ^iMu + hfj_^, we find, to linear order, 

8h 

= V2 V + 2d^t.v.) + 0{h'), (2.3) 

where 

Vu = -duhxx - dxhxy. (2.4) 

We see that the perturbation evolves according to the heat equation, accompanied by a 
small diffeomorphism.^ Just as with the scalar heat equation, we are restricted to Euclidean 
signature, since on a Lorentzian background the flow is infinitely unstable to modes with 
timelike momentum, and therefore ill defined. 

How unique is Ricci fiow? Suppose that we wish to define a flow equation for the metric 
that is local and diffeomorphism-invariant and that contains at most two second derivatives 



acting on the metric. There are only two terms we could add to the right-hand side of (2.1): 



-2R^^ + aRgf,^ + 2Ag^^. (2.6) 



Because it contains no derivatives, the A term does not change the behavior of short-wavelength 
perturbations, and therefore does not affect the locally diffusive nature of the flow. While it is 
useful when one is interested in flows that start or end on Einstein metrics with non-zero cos- 
mological constant, in this paper we will set A = for simplicity. On the other hand, having 
a non-zero a changes the right-hand side of ( |2.3| ), adding the terms a{d\dph\p — y'^h\\)5p_y, 
which has the effect of making different polarizations diffuse at different rates. A given 
perturbation can be decomposed into a part obeying dxdphxp — V^/iaa = 0, which evolves 
independently of a, and a Weyl part (proportional to 5pu), which evolves with a diffusion 
constant 1 — a{D — 1). In particular, unless a < 1/{D — 1) the flow is ill defined. In this sense 
Ricci flow is singled out; as we will discuss below, the value a = is also physically selected 
in the context of RG flow. It would be interesting to understand how general our results are 
for the other values of a for which the flow is well defined. 



""^For pure gauge perturbations, which are of the form /i^i, — for some vector field w^, the dif- 

feomorphism term cancels the Laplacian term, since the Ricci tensor remains zero. For some purposes (for 
example in order to rigorously prove the short-time existence of solutions, or for improved numerical stability 
in simulating the flow) it is useful to have a strictly parabolic flow equation, in other words one where all 
short- wavelength modes decay. An elegant way to do this, due to DeTurck po| and closely related to the 
harmonic gauge condition in General Relativity, is to fix an arbitrary connection F^^, on M (thereby breaking 
the diffeomorphism invariance), and define the vector field 

r = /^(rL-f^.), (2.5) 

where F^^^ is the connection compatible with gfj,^. One then adds 2V(^5y) to the flow equation, which only 
changes the flow by a A-dependent diffeomorphism. Expanding 2V(^^y) about flat space yields precisely 
— 2c)(^Vy), so that all short-wavelength perturbations evolve by the heat equation. As explained in Appendix 
A, a similar trick was employed in the numerical simulations of Ricci flow described in this paper. 
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Now let us review the basic global properties of Ricci flow. The fixed points are obviously 



Ricci-flat metrics. The generalization of the linearization (2.3) to the case g^y = g^i, + /i^j^, 
where g^i, is Ricci-flat, is 

= -Ai^h,, + 2V(^^,) + 0{h'), (2.7) 
where v,y = ^d,yh^x — V \h^y^ and Al is the Lichnerowicz operator: 

Al V = -V^Va V - 2i?M"/^«A (2.8) 

(everything is calculated with respect to the background metric g^u)- The linear stability of 
the fixed point is thus determined by the existence of negative eigenvalues of Al. (Al is non- 
negative on pure gauge modes, so any negative mode must be physical.) Al acts separately on 
Weyl {hfj_y{x) = e{x)g^y{x)) and traceless {h^^{x) = 0) modes. For Weyl modes it reduces to 
minus the scalar Laplacian, which is non-negative on a compact manifold. On the other hand, 
for traceless modes the Riemann tensor acts as a "potential" on the manifold (in the sense 
of the Schrodinger equation), with regions of positive curvature corresponding to negative 
potential energy. The operator may therefore have a finite number of "bound states", i.e. 
modes with negative eigenvalue. 

Another important property of Ricci flow, obvious from the definition, is that regions 
(and, within a region, directions) with positive Ricci curvature tend to collapse. This can 
lead to the formation of curvature singularities in finite flow time. The simplest example is 
to start with the round metric on the D-sphere. We will see a more complicated example in 
the flows we study in this paper. 

2.2 Boundaries 

In Section 4 we will consider Ricci flow on a manifold with boundary. We should therefore 
consider what boundary conditions are required to render the flow well defined. Let us locally 
divide the coordinates into ones that are parallel (x*) and normal (x"') to the boundary. In 
general for a diffusive flow one should impose a single boundary condition (such as Dirichlet or 
Neumann) on any variable / which is acted on by a second derivative normal to the boundary, 
i.e. where d^f appears, while for any variable acted on by fewer than two normal derivatives, 
no boundary condition should be imposed. In the Ricci tensor only the metric components 
along the boundary, gij, appear with two normal derivatives, so it is only on these that we 
should impose boundary conditions.^ For example, in the Dirichlet case we would fix the 



^ As explained in the previous subsection, the DeTurck flow is equivalent to Ricci flow up to a diffeomorphism 
but is strictly parabolic, so that all metric components appear with second derivatives in all directions. For 
this flow we will therefore need to impose boundary conditions on all components of the metric. The extra D 



boundary conditions (on gin and (?nn) are provided by the condition that (as defined in (2.5)) vanishes on 
the boundary. This is necessary for DeTurck flow to be equivalent to Ricci flow with the prescribed boundary 
conditions, which requires that the A-dependent diffeomorphism that relates them be the identity on the 
boundary. 
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induced metric 7jj on the boundary, while in the Neumann case we would fix its extrinsic 
curvature Kij (the latter case is treated in the paper These are the same kinds of 

boundary conditions one can impose on the hyperbolic time evolution of general relativity. 

In what follows, we will adopt Dirichlet boundary conditions. 7i is therefore defined as 
the space of Euclidean metrics on M with a fixed induced metric ^ij on dM. An important 
technical point is that, with these boundary conditions, we cannot make the usual decom- 
position of perturbations into Weyl and traceless parts. Specifically, a perturbation whose 
component doesn't vanish at the boundary cannot be expressed as the sum of an allowed 
Weyl and an allowed traceless perturbation. 



2.3 Ricci flow as gradient flow 

The fixed points of Ricci flow, namely the Ricci-flat metrics, coincide for D > 2 with the 
critical points of the Euclidean Einstein-Hilbert action, 

(5*0 is an arbitrary constant). This suggests that Ricci flow might be gradient flow of S. We 
will now show that this is the case. 

Recall that the general gradient fiow equation, on a space with coordinates and "en- 
ergy" function S, is 

The inverse metric G^^ is necessary to raise the index on the gradient. In our case is 
9fiu{x), the index A standing for both the point x in M and the component /iz^, and S is the 
Einstein-Hilbert action. But what should we take for Gab'^ The most general metric on Ti 
that is local and diffeomorphism-invariant (on M) and contains no derivatives is of the form 

G^dg^dg"" = j^^ ^ {dg>^udg\ + a{dg^^ ^f) , (2.11) 

where dg^^ = g^^dgxy. The only free parameter is a (except the trivial overall factor which 
has been fixed for convenience). With this metric we find 



(99 

-Gta^-d^ = -'^R,r^+<^R9,u, (2.12) 

where 



2q + 1 

a = — — — . (2.13 

Da + 1 



Gradient fiow thus establishes a one-to-one map between a and the parameter a from ( |2.( 
In particular, Ricci fiow (a = 0) is obtained when a = —1/2. We will simply refer to this 
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metric as Gab'-^ 

GABdg^dg^ = j^^ V9 [dg^dg^ - ^(d<7%)') • (2-14) 

Notice that Gab is not positive definite, as Weyl perturbations have negative norm. 
This is somewhat unconventional for a gradient flow, and implies that S does not necessarily 
decrease along the flow (hence we should avoid the terminology "gradient descent" ) . The fact 
that Gab is negative on Weyl modes is actually crucial for obtaining a well-defined gradient 
flow from the Euclidean Einstein-Hilbert action. (In fact, notice that more generally the 
condition a < 1/{D — 1) mentioned below ( |2.6| ) requires a < —1/D, i.e. requires the Weyl 
modes to have negative norm.) The reason is that short- wavelength Weyl modes decrease the 
value of S. A conventional gradient descent would be therefore be unstable to all such modes. 
However, if they have negative norm then the flow runs uphill in those directions, rendering 
it well behaved. 

The Weyl modes' negative norm plays a similar role in Euclidean Quantum Gravity, 
where one integrates over Ti. In that context, the fact that these modes decrease the 
action, which is therefore not bounded below, threatens to make the path integral ill defined. 
The standard cure is to rotate the contour of integration for each Weyl mode to lie parallel to 



the imaginary axis |23|. In analogy with the path integral for a particle or string moving in a 
Lorentzian background, where the contours for the timelike target space direction are similarly 
rotated, this prescription is quite natural if one takes Gab as the metric on TC. We know 
of no argument for Euclidean Quantum Gravity that requires this metric to have a = —1/2. 
This value naturally gives rise to the de Bonder gauge condition when removing pure gauge 
modes from the one-loop determinant, as we will review in subsection 3.1. However, it would 
be interesting to study the path integral and its gradient flows arising from different choices 
of a. 

2.4 Relation to sigma-model RG flow 

A two-dimensional sigma model is defined, at a given energy scale A, by a target space 
manifold M equipped with a metric g^u{x), a two-form i3^,y(x), and a dilaton $(x). Let us 
set Bfj^i, = (which is equivalent to imposing parity symmetry). To first order in the coupling 
a' , the running of the metric and dilaton is given by 



^ = - 4V^9,$, II = + - 2(5^$)^ (2.15) 



where A = — ^a'lnA. (We have written the equations for the bosonic sigma model, but the 
only change for the super symmetric versions at this order in a' is the coefficient of 1/a' in the 



^In a previous version of this paper, we mistakenly referred to this metric as the "DeWitt" metric. In fact, 
in the DeWitt metric, which plays an important role in the Hamiltonian formulation of General Relativity 
p^, a takes the value —1. We would are very grateful to S. HartnoU for pointing this out, as well as for several 
other very useful comments on this subsection of the paper. 
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dilaton equation.) The flow equation for g^u is of the form ( p. 2]) , with = — 2(9^<I>. Therefore, 
irrespective of how the dilaton evolves, the metric evolves according to Ricci flow, up to a 
A-dependent diffeomorphism. (However, one as to be careful: because of the dilaton a fixed 
point of RG flow is not necessarily a fixed point of Ricci flow, and vice versa.) Alternatively, 
one can consistently set the dilaton to a constant, in which case g^y evolves strictly according 
to Ricci flow. 

Just as Ricci flow is gradient flow of the Einstein-Hilbert action, the RG flow defined by 
( |2.15| ) is gradient flow of the spacetime effective action of string theory, 

S[g,,{x), $(x)] = I Vae-'" (-^ + R + <d,^f) - \ ( V^e-'^K, (2.16) 

Jm \ 6a J 1^ JdM 

with respect to the metric 

4^ / [d9''ud9\ - \{dg^, - M^?^ . (2.17) 

Note that the timelike direction is now a particular mixture of the dilaton and the Weyl mode 
of the metric. 

The RG flow of two-dimensional quantum field theories is a subject of physical interest 
in its own right. Recently it has also become a popular tool for exploring the off-shell con- 
figuration space of string theory, for example in the study of tachyon condensation.^ We will 
not include a dilaton in the analysis that follows, but since as explained above the dilaton 
doesn't enter into the running of the metric, our work can be seen as an off-shell exploration 
either of pure gravity or of string theory. 

3. Gravity in a cavity 

In this paper we will study Ricci flow on 4-manifolds with boundary x S^, employing 
Dirichlet boundary conditions with the boundary metric fixed to be a round sphere of radius 
R times a circle of circumference (3. These boundary conditions arise in two distinct physical 
contexts. The first is the canonical ensemble for 4-dimensional gravity in a spherical box, 
with the circle representing the imaginary time direction. The second is the superspace for 
5-dimensional gravity Kaluza-Klein reduced on a circle, again in a spherical box. 

We will begin by reviewing the canonical ensemble for gravity in a box. We will then 
explain how the mathematical features of this system are to be physically reinterpreted in 
the Kaluza-Klein context. 

3.1 Thermal interpretation 

The canonical ensemble for gravity in asymptotically flat space at finite temperature is ill de- 
fined due to infrared instabilities, both perturbative (Jeans instability) and non-perturbative 

*See and references therein for a discussion of the relationship between world-sheet RG flow and 
spacetime dynamics of closed string tachyons. 
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(nucleation of black holes). An infrared regulator is therefore necessary, the simplest one 
being the spherical box we employ here. Other regulators are possible; for example one can 
introduce a negative cosmological constant and work in asymptotically anti-de Sitter spaces. 
The thermodynamics of the system doesn't depend qualitatively on which regulator is used. 
What is described in this subsection was largely the work of Gross, Perry, and Yaffe |]l|. 
Hawking and Page |^5[, and York [^], except for the calculations of the negative mode for 
the small black hole (resulting in figure ^) which are new. 

We consider pure gravity in a rigid spherical box of radius R, which is kept in contact 
with a heat bath at a temperature Note that it is the boundary that is kept at that 

temperature; even though the system is in thermal equilibrium, observers at different points 
inside the box may experience different temperatures due to gravitational red- and blue- 
shifting. According to the prescription of Euclidean Quantum Gravity, the partition function 
for this system is obtained by integrating e"'^'^' over all Riemannian manifolds with boundary 
metric Sj^ x S]^. This involves both summing over topologies and integrating over metrics.^ 

The path integral is usually treated by the saddle point method. The number of saddle 
points (i.e. Ricci-flat metrics) depends on the dimensionless ratio b = P/{2ttR). For any 
value of b, we have hot flat space, Sj^ x Bj^, where Bj^ is a flat 3-ball of radius R. For 
b < 6crit = 4/3\/3 « 0.77 (in other words, at sufficiently high temperature), there are two other 
solutions. Both of them have topology B^ x S'^, and both have the Euclidean Schwarzschild 
metric, but with different values of the horizon radius tq: 

ds^ = Arl (l - ^) dr^ + (l - ^)~' dr^ + r^d^l (3.1) 

The radial coordinate r goes from tq (at the horizon) to R (at the boundary). The imaginary 
time coordinate has periodicity 27r, and the factor of Ar^ in the metric ensures smoothness at 
the horizon. For the circle to have the right length at the boundary, we must have 

/? = 47rroyi-^, (3.2) 

which fixes ro in terms of (5 and R. For b < 6crit there are two solutions, with rg greater 
than and less than |i? respectively, which we refer to as the small and large black holes. 
Although it has not been proven, it is believed that these are the only saddle points for these 
boundary conditions [^], and this is what we will assume. In the limit that we remove the 
cutoff (R oo, (3 fixed), the small black hole goes over to the usual asymptotically flat 
Schwarzschild black hole (ro — > /3/47r), while the large black hole disappears (ro oo). Hot 
flat space persists as well in this limit, of course. 

The system's partition function, and therefore its thermodynamical properties, are dom- 
inated by the saddle point with the lowest value of the action (^]^). Choosing the constant 

^In principle this includes disconnected manifolds. However, in this case the given boundary is connected, 
so any manifold satisfying the boundary conditions can be divided into the connected component containing 
the boundary and the other components. The other components factor out of the path integral, multiplying 
the partition function by a constant factor. Therefore it is sufficient to consider connected manifolds. 
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Figure 1: Action versus inverse temperature parameter h = f3/{27TR) for tlic three saddle points. For 
b > fccrit only hot flat space is allowed, while for b < 6crit there are also small and large black hole 
solutions. At 6hp there is a first-order phase transition, the Hawking-Page transition: for b > 6hp 
flat space has the lowest action and therefore dominates thermodynamically, while for b < 6hp the 
large black hole dominates. The dashed lines indicate the three values of b which will be used in later 
figures to display the qualitative behavior of the simulated Ricci flows. For this and all other flgures 
we have chosen units such that R = 1. 



So to conform to the convention that flat space has zero action, the Schwarzschild solutions 
have action 

The actions of the three saddle points are plotted as a function of b in figure || (for this and 
all other figures we choose units such that R = 1). We see that, for all temperatures where it 
exists, the small black hole has the largest action of all the saddle points. At low temperatures, 
specifically for b > 6hp = 16/27 ~ 0.59, hot flat space dominates, while at higher temperatures 
the large black hole dominates. The first-order phase transition separating these two regimes 
is known as the Hawking-Page transition, by analogy with the same phenomenon in AdS |25]. 

So far we have neglected the one-loop determinant in the saddle-point approximation to 
the path integral (not to mention the higher-loop corrections, which require an ultraviolet 
cutoff, such as string theory, to define). This correction is indeed negligible at small h (i.e. 
when the system is much larger than the Planck length), as long as all the eigenvalues are 
positive; otherwise the saddle point approximation doesn't make sense. The operator in 
question, taking into account the rotation of the contours for the conformal modes, is 
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Figure 2: Lowest eigenvalue [i of the Lichnerowicz operator Al on the black hole background (3.1). 
The small and large black holes are separated by the value r^jR = |, where the eigenvalue passes 
from negative to positive. At this point the mode becomes tangent to this one-parameter family of 
solutions. 



Acting on a perturbation h^p^ yields precisely (minus) the right-hand side of (2.7) 

(since the latter is just dc{—G^^dBS)h'^, and of course ObS = at the saddle point). This 
operator has infinitely many zero eigenvectors corresponding to pure gauge fluctuations, of 
the form hfi^^^ = V(^t<;^), which must be removed before the determinant is calculated. 
Since is self-adjoint (with respect to Gab) the physical eigenvectors are orthogonal 

to the pure gauge ones, and therefore correspond to fluctuations obeying the de Bonder 
gauge condition Vf.h'^u - \duh = (since GABh^s,ngeh^ = I \/9^ (,fiWy){h^'' - ^hgf"") = 
— f ^/gWl,{S/ ^h^^ — ^d'^h)). Restricting to the orthogonal complement of the gauge zero 
eigenvectors leaves these physical fluctuations on which A^c is just the Lichnerowicz operator 
Al. (Equivalently, one can follow the Faddeev- Popov procedure and add a gauge-fixing term 
and ghosts 

On the hot flat space background, Al is positive since the Riemann tensor vanishes and 
we have imposed Dirichlet boundary conditions.^ What about on the black hole backgrounds? 
In figure |2|, the lowest eigenvalue of Al is plotted against tq. For the large black hole (ro > f-R) 
it is positive, so like flat space it is indeed a good candidate for the thermal ground state of 
the system. In the regime where each has the lowest action, that action (divided by /?) is the 
classical approximation to the system's free energy. The small black hole, on the other hand, 
has a negative eigenvalue mode, the cavity generalization of the Gross-Perry- Yaffe instability 
of an asymptotically flat black hole [||. As in the asymptotically flat case, this negative mode 
is mainly localized near the horizon, and preserves the U{1) x SO{3) isometrics of the black 

®At extremely high temperatures there are more subtle instabilities which we do not consider here. When 
(3 is smaller than about gI^'^'R^^^ the gas of thermal gravitons suffers from the Jeans instability Within 
string theory, when f3 is smaller than the string scale there will be an unstable winding state, indicating a 
Hagedorn phase transition [E8|. 



hole (i.e. only depends on the radial coordinate). In past literature there is some confusion 
over the calculation of this mode in a cavity [p9| , p6| , The important point is that, unlike in 
the asymptotically flat case, the mode is not traceless, since the cavity boundary conditions 
(constant Euclidean circle and 2-sphere radius) explicitly couple the trace and traceless modes. 
When this is taken into account one finds the curve in figure ^. (See Appendix A for details of 
the calculation.) In particular, as expected, one finds that the instability disappears exactly at 
the transition from the small to the large black hole, and therefore at the transition from local 
thermodynamic instability to stability. As pointed out by Hawking and Page ||2^, the zero 
mode that appears corresponds to slightly changing the black hole's radius without altering 
the boundary condition. 

What is the meaning of this negative mode? It has been argued to be related to the fact 
that the corresponding Lorentzian solution is locally thermodynamically unstable since it has 
a negative specific heat [||, ^. The negative mode also means that the Euclidean solution 
can serve as a bounce, representing thermal fluctuations that take hot flat space to the large 
black hole and vice versa, allowing the system to find its true ground state whichever side of 
the Hawking-Page transition it's on [Q, |2^. So although it is never the lowest-action saddle 
point, the small black hole nevertheless plays an important physical role in the system's 
dynamics. 

3.2 Kaluza-Klein interpretation 

The three 4-dimensional Euclidean Ricci-flat geometries discussed in the previous subsection 
have an entirely different physical meaning if we interpret the circle direction, which was 
the imaginary time in the thermal context, instead as a Kaluza-Klein circle, and add a real 
time direction. The resulting 5-dimensional Lorentzian geometries are static solutions to the 
5-dimensional Einstein equation with boundary R x 5"^ x S]^. Note that, whereas in the 
4-dimensional context these solutions are believed to be the only ones, in the 5-dimensional 
context there are certainly other static solutions, such as the uniform black string, which is 
(4-dimensional Lorentzian Schwarzschild) xS"^. This is because there is an additional degree 
of freedom present, namely the time-time metric component. 

What do these 5-dimensional solutions look like? First note that the ADM energy of 
each is equal to the Euclidean action of the corresponding 4-dimensional solution, plotted 
in figure |l[ Furthermore, its linear stability is controlled by the 4-dimensional Lichnerowicz 
operator, since each mode of that operator can be put on shell by dressing it with a time 
dependence e^*v^*, where n is the Al eigenvalue. Let us fix /?. For all values of R we have 
"hot flat space" , which in this context is the Kaluza-Klein vacuum. As pointed out by Witten, 
although it is perturbatively stable, in the limit i? — > oo the Kaluza-Klein vacuum is unstable 
due to nucleation of bubbles of nothing |3C]. At large but finite R we might therefore expect 



to find a ground state with negative energy; indeed we find the "large black hole", which in 
this context is a static bubble of nothing that almost fills the container, but leaves a shell 
of "something" around the edge. This is also perturbatively stable. It seems likely that, for 
b > 6hp, the true ground state is the Kaluza-Klein vacuum, while for b < 6hp it is the large 
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bubble of nothing. Finally, we have the "small black hole" , or small bubble of nothing, which 
has higher energy than either the Kaluza-Klein vacuum or the large bubble, and is unstable/ 
In the context of string theory, such a gravitational instability is classified as a closed 
string tachyon. It is always interesting to study the world-sheet RG flow seeded by the 
vertex operator of such a tachyon, and to compare it to the system's dynamical evolution. 
The classical evolution of the bubble perturbed by its unstable mode has in fact already 
been studied by Sarbach and Lehner using numerical simulation ||^, |8|. Those simulations 
were performed without the box, which simplifies the analysis because it allows radiation to 
escape to infinity and the energy thereby to decrease. It was found that, when perturbed 
in one direction, the bubble grew without bound, while in the other it settled down to the 
black string. Apparently cosmic censorship prevents the system from classically reaching its 
local ground state, the Kaluza-Klein vacuum, because that would entail the appearance of 
a naked topology-changing singularity. (Of course quantum mechanically the black string 
would presumably evaporate, eventually settling down to the vacuum.) We will see in the 
next section that RG fiow (which is well approximated by Ricci flow when the system is much 
larger than the string scale) does not suffer from such prudishness. 

4. Where does the small black hole flow? 

Let us call the negative mode of the small black hole h'^^^ . If we perturb the small black hole 
by this mode and evolve the metric by Ricci fiow, then, as long as it remains sufficiently small, 
the mode will grow exponentially. How does the metric evolve once it leaves the linear regime? 
At late times does it converge to one of the other saddle points? These are the questions we 
will answer in this section. Note that there are two different fiows to follow, depending on 
whether h'jj^^ is added with a positive or negative coefficient. These are physically distinct; 
in one direction the mode increases the horizon size and in the other decreases it. We will 
fix the convention that h^^^ increases the horizon size. We performed numerical simulations 
of both of these Ricci fiows at several values of the parameter 6, below, at, and above the 
Hawking-Page value 6hp- Details of the simulations are given in Appendix A. In this section 
we summarize the main results. 

First, however, we would like to point out two obvious generalizations of our work. The 
first is to dimensions D > 4. The second is to the small black hole in anti-de Sitter space, in 
which case the fiow equation must be modified to include the cosmological constant term: 



In both cases the system is so closely analogous to the one we studied that it seems reasonable 
to guess that the behavior of the fiows would be qualitatively the same. 

We begin with the fiow seeded by the perturbation +h^y^ ■ Since the perturbation does 
not break the black hole's U{\) x SO{2>) isometry group, the geometry enjoys this symmetry 



dX 



2R^^ + 2Ag, 




These solutions are the analogs for zero cosmological constant of the asymptotically AdS solutions found 
by Copsey and Horowitz jsij. 
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Figure 3: Euclidean time radius T (red) and sphere radius S (purple) against proper radial coordinate 
p, for the flow of the small black hole perturbed by at & = 6hp ~ 0.59. The boundary is at 

p — Q. Snapshots are drawn at intervals of A of 0.05. The heavy blue curves show those functions 
for the unperturbed small black hole. The heavy black curves show them for the corresponding large 
black hole, to which the flow clearly asymptotes at late times. 

throughout the flow. In principle we would consider an infinitesimal perturbation. In practice 
for the initial data we perturbed the small black hole such that the horizon was increased in 
radius by one percent. This was sufficiently small that for some time the flow remained well 
described by the linear theory. 

To characterize the geometry it is sufficient to know the radii of the time circle and sphere, 
as functions of proper radial distance; in other words, to eliminate the gauge ambiguity the 
metric may be put into the form 

ds^ = T{pfdT'^ + dp^ + S{pfdal, (4.2) 

setting yO = at the cavity wall (and p < in the cavity interior). T{p) and S{p) are plotted 
in figure ^ for a representative value of 6, namely 6hp- We see that the metric tends at large 
A to the large black hole solution. Let pn to be the value of p at the horizon. In figure ^ the 
horizon radius, i.e. ^(ph), is plotted as a function of A for three different values of b. Again, 
in each case it asymptotes to the value for the large black hole. The initial growth of the 
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Figure 4: Horizon radius, S'(/3h), against A for small black hole perturbed by +h'j^^^ , a.t b ~ 
0.40,0.59,0.70 (the middle value equals 6hp)- The horizontal (red) dashed lines show the horizon 
radius for the corresponding large black hole, which the metric asymptotically converges to in each 
case. The curved (orange) dashed lines give the initial exponential growth of the negative mode. 

small perturbation is governed by the linear theory, growing exponentially as e~^^, with /u 
given in figure The same behavior was found for all values of b tested: the metric relaxes 
to the large black hole solution. 

4.1 Singularity and surgery 

The behavior was quite different for the flows seeded by the perturbation —h^^^- For every 
value of h tested, after a finite flow time the horizon collapsed to a point, giving rise to a 
curvature singularity and stopping the flow. An example is shown in the top panel of figure 
|. In fi gure y, the horizon radius is plotted against A for the same three values of b as shown 
in figure ^. More details of the local model for the singularity are given in appendix B. 

The spontaneous formation of singularities is generic in Ricci flow, and has been much 
studied by mathematicians, particularly in the case D = 3. A common type of singularity that 
arises there, and is analogous to what we observe here, involves an S*^ fibered with varying 
radius over an interval (in a region of the manifold). The fibers' positive curvature tends to 
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Figure 5: (Top) T{p) (red) and S{p) (purple) for the flow of the small black hole perturbed by 
—h^^^: at 6 = 6hp. Snapshots are drawn at intervals of A of 0.01. The metric flows to a singularity, 
where the horizon shrinks to zero size. (Bottom) The metric after surgery, and the continuation of 
the Ricci flow. The metric asymptotically tends to flat space. 

make them shrink, and typically in finite time one of them will collapse to zero size, giving rise 
to a curvature singularity. In order to be able to continue the flow, mathematicians perform 
"surgery" on the singular manifold: a small neighborhood of the singular point is excised and 
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Figure 6: Horizon radius, S{ph), as a function of A for the flow of the small black hole perturbed by 
—h^^^- In all such flows the horizon collapses to zero size in finite time. The dashed lines show the 
initial exponential growth of the negative mode given by perturbation theory. 

replaced with a smooth metric. To avoid having the singularity re-appear, the topology is 
changed. Specifically, a piece of the base is removed, and smooth caps are sewn onto the two 
ends, as shown in figure |^. As a result of this process, the manifold (or at least this part of 
it) becomes disconnected, and the S"^ fibers, which were previously non-contractible, become 
contractible. 

The singularity that appears in the Ricci flow of the black hole is quite analogous. Again, 
we have an S"^ fibered over a base, which in this case is a disk. And again, the fiber over a 
point on the base collapses to zero size after a finite flow time. We are therefore inspired by 
the above example to perform a similar surgery in order to continue the flow. Specifically, 
we excise a small disk from the base in a neighborhood of the singular fiber. The base thus 
becomes an annulus, and on its new inner boundary the fiber radius is made to go to zero 
in order to produce a smooth total space with topology x B^; the 5^, which was non- 
contractible, is now contractible, and vice- versa for the S^. The change to the metric is 
localized to a small distance e from the singularity, and is required to respect the geometries 
C/(l) X SO{'i) isometry group. A cartoon of it is shown in figure |8|. The resulting geometry 
can be seen in the example of figure |5| in the first snapshot of the bottom panel. Details of 
how the surgery was carried out are given in Appendix B. 

Those who are interested in Ricci fiow chiefly insofar as it approximates the RG flow 
of sigma models might ask the following question: What happens to the RG flow when the 
Ricci flow hits the singularity? When the curvatures become large in string units, the a' 
corrections to the beta functions ( |2.15| ) become important, and Ricci flow is no longer a good 
approximation to RG flow. However, in general this will not be sufficient to stop a singularity 
from forming. This is shown by a simple example, namely the round metric on (otherwise 
known as the 0{N) model, with N = D + 1). On the Ricci fiow side, the sphere collapses to 
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Figure 7: An S'^ fibered over an interval will typically collapse under Ricci flow, producing a curvature 
singularity. This can be cured, allowing the flow to continue, by a local topology change known as a 
surgery. 

zero size in finite time. On the RG flow side, the model confines at a finite energy scale, and 
is massive below that scale. In the case of a sphere fibered over a base, it seems reasonable to 
guess that, similarly, when the fiber over some point on the base collapses ("confines"), it is 
removed from the target space, whose topology would change accordingly.^ In other words, 
we are proposing that, when the Ricci flow hits the singularity, the RG flow dynamically 
performs the surgery we have described above. The fact that this surgery is the unique local 
topology change that respects the manifold's isometrics means that the RG flow could hardly 
do otherwise. 

After the surgery, the manifold has the same topology as hot flat space, and under Ricci 
flow the metric relaxes to that saddle point without encountering any further singularities. 
This can be seen in the bottom panel of figure |; the same behavior was observed for the 
other values of h tested. A key assumption is that, provided the metric deformation in the 

*This is similar to the picture of |^^, but in that case the fiber is an and the confinement is driven by 
the condensation of a winding tachyon. 
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Figure 8: When perturbed by —h'j^^^, the smah black hole geometry (which is a 2-sphere fibered 
over a disk) evolves under Ricci flow until the horizon (tire fiber over the center of the disk) collapses, 
giving rise to a singularity. The surgery we perform on the singular manifold involves excising a small 
disk from the base, making it an annulus, and making the sphere go to zero size on the inner boundary 
so that the resulting total space is smooth and topologically x B^. We do this in such a way that 
the U{1) X S0{3) isometrics of the geometry are preserved. 

surgery is restricted to a region of size e, for small enough e the details of the flow following 
the resolution will be independent of the precise nature of this deformation, and of the actual 
value of e. This is to be expected since the diffusive nature of Ricci flow quickly removes 
the fine details of the small resolved region. While we do not have a rigorous mathematical 
argument that this is the case, we found it to be true in our simulations, as shown in Appendix 
B. (The value of e used to generate the singularity resolution in figure ^ was 0.07). 

5. A novel free energy diagram 

In the last section we saw that Ricci flow defines paths in the space of metrics connecting 
the small black hole to the large black hole and to hot flat space, in the latter case via a 
topology change. Since Ricci flow is gradient flow of the Einstein-Hilbert action ( |2.9D , it is 
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Figure 9: Action S versus flow time for the flows from the small black hole to hot flat space (left) 
and the large black hole (right), at 6 = 0.40, 0.59, 0.70 (top, middle, bottom). The black lines indicate 
the topology x 5^, and the blue the topology x B^, their join in the left plots being the point of 
singularity formation where the action remains continuous (up to a small numerical error in the case 
b = 0.40). 

interesting to plot the action along the flow. This is done in figure ^ for the two flows at 
the three representative values of b used in previous figures. We would like to point out two 
interesting features of these plots. 
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Figure 10: S/(3 against proper distance along tlie one-parameter family in the space of metrics 
defined by Ricci flow. In the left-hand (blue) part of each curve the topology is x B^, while in the 
right-hand (black) part the topology is x S^. As explained in the text, for these curves the value 
of S/ (3 can be interpreted as a free energy. 

First, in the flow to hot flat space, the action stays finite when the singularity forms and 
is continuous across the topology change; although the Ricci scalar blows up on the horizon 
when it collapses, it remains integrable. The local singularity model described in Appendix 
B indicates that the divergence in the Ricci scalar is entirely due to the shrinking 2-sphere, 
so that i? ~ 1/5"^ at the horizon. This is compensated by the factor of the 2-sphere volume 
in the measure ^ ~ TS'^, and since T ^ at the horizon too, there is no contribution 
to the action from this singular point. From the point of view of the action, therefore, the 
singularity is harmless, and the surgery can be thought of as an infinitesimal change in the 
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metric. 

Second, the action is a monotonicahy decreasing function of A. This is not a straightfor- 
ward consequence of the fact that Ricci flow is gradient flow of the action, since the metric 
Gab is not positive definite (which is why we have used the term gradient flow, rather than 
gradient descent). The monotonicity of S for these particular flows shows that they move 
through the space of metrics always in a spacelike direction; in other words, at every point 
along the flow we have 

8S_^ Bg-^Og" _ 1 ( „„„ 1 . 



Hence the proper distance along the curves is well defined; it is given by 



(5,2) 

In figure p!o|, we plot S against q, putting the two flows starting at the small black hole back 
to back where they join smoothly. Note that the saddle points are a flnite distance away from 
each other along these curves. 

Figure |l^ bears a strong resemblance to a free energy diagram. Can it be interpreted 
as such in a physically meaningful way? To clarify this question, let us flrst briefly review 
how free energy diagrams are conventionally deflned. Suppose one is interested in knowing 
the expectation value of some macroscopic observable a (for example, the magnetization of a 
ferromagnet). To do this, one calculates the free energy /(a) in the sub-ensemble where the 
value of a is flxed; this calculation can, for example, be done in a saddle point approximation. 
The expectation value of o, and the actual free energy F of the system, are determined 
by integrating e~^^'^°'^ over o, which, again in a saddle point approximation, amounts to 
minimizing /(a). In other words, in calculating the partition function the degrees of freedom 
are integrated over in two steps: flrst everything but a, then a itself. The free energy diagram 
appears between these steps. 

Let us illustrate these ideas for the case at hand, gravity in a cavity; our discussion 



follows York |26]. We choose as our macroscopic observable a the horizon area. Working on 
the topology x 5^, the horizon area is flxed by adding a Lagrange multiplier term to the 
action: 

Sa[g] = S[g]+X{A[g]-a). (5.3) 

This term introduces into the Einstein equation a source term in the form of a two-dimensional 
delta function on the disk at the location of the horizon, leading to a conical singularity. The 
saddle point of the new action is therefore the Euclidean Schwarzschild metric ( |3.1| ), but 
with ro determined by A-nr^ = a rather than by ( [3.21) , and with the periodicity of r equal to 
(3/ {2rQ^l — tq/R) rather than 2it. The action of this solution is easily calculated, giving a 
free energy of the familiar form "energy minus temperature times entropy" : 



f = ^ = A 




Af3R 



(5.4) 
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Figure 11: Conventional free energy diagram for gravity in a cavity, in which the free energy / 
is calculated in sub-ensembles where the horizon area a is fixed. In the saddle point approximation, 
f = S/P, where S is the action of the Schwarzschild metric with horizon area a and a conical singularity 
in the transverse directions. The three curves are for the same values of b as in the previous figure. 



Following Whiting and York we plot / in figure 11, where the three saddle points are 
easily identified. Hot flat space is at o = 0; in fact the topology change from x 5^ to 
X B"^ occurs at that very point on the free energy diagram (which to some extent excuses 
the rather singular behavior of / there). 



Let us now return to our putative free energy plot, figure |10|, derived from the Ricci 
flows. S/P, which is plotted for each point g'^iX) along the flow, can justifiably be called 
a free energy if it is the saddle point approximation to the partition function of some sub- 
ensemble. We will now show that it is. First note that the gradient of S, restricted to the 
directions orthogonal to the curve g^{X), vanishes, since the tangent vector to the curve 
is in the gradient direction. Furthermore the second-derivative operator A'^c (defined in 
( |3.4D ) restricted to those directions is positive, since its only negative direction is along the 
curve. Therefore the saddle point approximation to the path integral along the directions 
orthogonal to the curve is given by the value of the action on the curve. This defines the 
sub-ensembles for which figure ^ is the appropriate free energy diagram. It is obviously 
crucial for this construction that g'^iX) is defined by a gradient flow of the action. In fact 
the same construction could be applied to any system equipped with a natural metric on the 
configuration space. The major difference between this construction and the conventional one 
is that this one does not require singling out a particular macroscopic observable. 

We note that York's prescription of fixing the horizon area is motivated by the Lorentzian 
dynamics of an evaporating black hole, since the horizon size is the slow variable in that 
process. From the point of view of the canonical Euclidean path integral, York's construction 
is less natural, as seen by the conical singularity that arises in the Euclidean geometries. 
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We conclude by commenting that in the context of the AdS-CFT duaUty it is inter- 
esting to consider the constrained free energy of asymptoticahy AdS spacetimes with and 
without horizons. On the field theory side of the correspondence, at weak coupling, one 
can naturally compute such an off-shell free energy [^^, which has a similar appearance to 
ours. At strong coupling, this is expected to be similar, with the extrema corresponding on 
the gravity side to exactly the small and large AdS-Schwarzschild black holes, and thermal 
AdS, and the Hawking-Page phase transition being interpreted in field theory as the con- 
finement/deconfinement transition. We note however, that the constraint in the field theory 
is placed on the expectation value of the trace of the Polyakov loop. On the gravity side, 
this presumably translates into a constraint on the area of a string world sheet wrapping the 
thermal Euclidean time circle at the asymptotic AdS boundary. Hence this is a different con- 
straint than the one imposed by York; it would be interesting to find the relevant constrained 
off-shell geometries.^ 
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A. Details of Ricci flow simulation 

In practice we choose units so that R = 1. Then the small black holes are given by ( |3.1| ) for 
ro < |, and the large black holes with | < tq < 1. Our boundary conditions on the box 
wall are fixing the sphere radius {R = 1), and fixing the Euclidean time period P = 27r6. We 
start with the small black hole metric, perturbed by a suitably small h*^^^ which satisfies the 
above boundary conditions. For the data presented here, the amplitude of the perturbation 
was chosen to change the horizon radius by one percent. The initial perturbation it/i^^^ 
respects the U{1) x S0{3) isometry of the black hole, and hence the entire flow will share 
this isometry. We exploit this, using the metric, 

ds' = e^^'b^dr^ + e^^dr^ + e^^d^l (A.l) 

where r has period 27r. Since we initially are interested in the topology x 5^ then A,B,C 
are finite functions of r and the flow parameter A, and we take < r < 1. At r = 1 we take 
A = 0,C = appropriate to our finite box boundary condition. The new 'boundary' at r = 




for related work on this system. 
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is fictitious, since we know the Euclidean black hole geometries are smooth manifolds without 
an interior boundary. Regularity of the Ricci flow equations implies A^^ = B ^ = C,r = 
at this origin. In particular, one finds as expected that {A/C)\r=o is preserved under the 
flow, and hence an initially smooth 4-d solution remains smooth, without developing conical 
singularities. 

For asymptotically flat spacetime we are free to make any choice of .^^ in ( |2.2D . However, 
on a manifold with boundary we must take care. We have covered our manifold with coor- 
dinates such that the boundary is at r = 1. Hence Ricci flow with = is only equivalent 
to flow with non- vanishing provided .^^ = at r = 1. Otherwise the coordinate location of 
the boundary will change during the flow. 

Subject to = at the box wall, we are free to specify any in ( |2.2| ). We make the 
choice, 

e = -Asx^ + Asx^U,ij,c=o (A.2) 

where are the coordinates and As is the scalar Laplacian. This choice is similar to 
the DeTurck flow, and ensures strong parabolicity — in analogy with the harmonic coordinate 
choice which yields strongly elliptic Euclidean Einstein equations. The flow equation becomes, 

-^g^^ = As5m^ + F^ia, V(7) (A.3) 

where on the right-hand side we have only explicitly displayed the second derivative terms, 
which are simply given by the scalar Laplacian. This parabolic flow requires boundary con- 
ditions for all components of the metric, namely A,B,C, at the cavity wall. There we take 
A = 0, C = 0, since we fix the induced metric there. The only non-zero component of is, 

= -drA + drB - 2drC (A.4) 

and hence we see the boundary condition that ^j. = at r = 1 translates into a boundary 
condition of the gradient of B. Hence A,B,C all have boundary conditions on their value, 
or gradient, compatible with the parabolic fiow. 

It is interesting to note that the 4-d Schwarzschild metric takes the elegant analytic form, 



(A.5) 



in the coordinate system where = everywhere, and furthermore the derived metric 
functions A,B,C are finite as we require. Hence we take this as a background metric. We 
compute using the Lichnerowicz operator in this background, also in the gauge = 0. 

As usual, the resulting coupled ordinary differential equations are best solved as a shooting 
problem. Note that the condition = together with A = C = implies that we cannot 
expect B = at the cavity wall. As mentioned in the main text, this means the mode cannot 
be traceless there, which clearly would also require B = 0. Hence the cavity wall couples the 
trace and traceless perturbations, unlike in the asymptotically flat or AdS computations. 
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Having computed this, we use it to perturb the small black hole with the amplitude 
chosen as described above. This forms our initial data for the Ricci flow, and ensures our 
condition = is satisfied at the box wall. 

We simulated the full Ricci flows for 

b = 0.10, 0.20, 0.30, 0.40, 0.59, 0.62, 0.64, 0.66, 0.68, 0.70, 0.72, 0.74. (A.6) 

The value b = 16/27 ~ 0.59 gives the middle dotted line in figure || where the free energies 
of hot flat space and the large black hole are equal. The values 0.40 and 0.70 give the other 
dotted lines in this figure. Properties of the flows for these three values are displayed in the 
various figures here, but we obtained qualitatively similar behavior for all values of b above. 

We used a uniform lattice discretization, and solve the Ricci flow using the implicit 
Crank-Nicholson method since it has diffusive character. We performed flows at various 
resolutions and confirmed the correct convergence to the continuum. Data presented here use 
dr = 0.0025, i.e. 400 points. 

After singularity resolution for the —h*^^^ flows, we are interested in the topology x 
where we use the metric, 

ds' = e^h^dr^ + e'^'dr^ + e^^r^dO^ (A.7) 

and take finite A,B,C, again with < r < 1. Now regularity of the Ricci flow equations 
at r = requires the boundary conditions B = C and A^r = B ^ = C^r = at the origin, 
ensuring again a regular 4-d Euclidean geometry. As for ( [A.l| ) we make the choice, 

e = -Asx'' + ^sx>'\A,B,c=o^ (A.8) 

again ensuring a strongly parabolic flow. As above, we require j4 = C = at the cavity wall, 
and obtain a boundary condition for B from the condition that vanishes there. 

B. More on the singularity and surgery 

In figure |l^ we plot the curvature invariant {Rap^vR°'^^'^Y^^ for the flow generated by —h^y^- 
We see that, as expected, the shrinking of the 2-sphere at the horizon leads to an increasingly 
singular curvature there. This nature of the singularity that forms may be seen by plotting 
T, S in the vicinity of the horizon, as in figure p. We see that since the geometry is smooth 
before the singular flow time, the gradient of T is fixed at the horizon. The interesting part 
of the geometry is then the shrinking 2-sphere, and we see that whilst the size of this at 
the horizon decreases towards the singularity, the rate of expansion of the sphere away from 
the horizon appears to remain zero as the singularity is approached. Hence the local model 
for this singularity is simply x S*^. Note that the only other option we know of for the 
singularity model would have been a cone with base x S"^. 

We wish to perform surgery on the geometry just before the singularity is reached. We 
do this by taking the geometry at some flow time Aq before the singularity, with functions 
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Figure 13: T{p) and S{p) plotted against proper distance from the liorizon, {p — pn), in tlie vicinity 
of the horizon as the singularity is approached. Snapshots are taken at intervals of A of 0.002. 



A, B, C in A.l, and moving to the topology |A.7 , with, 



e^^ = 
e^^ = 

e^^ = 


(e2 + r2)e2^ 


(B.l) 
(B.2) 

(B.3) 






e = 


e \r=0- 


(B.4) 



with, 



Note the resulting geometry with topology x is smooth, and for r ^ e the geometry 
before and after surgery is the same, the modification only occurring for r ~ e. The size 
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Figure 14: Comparison of two different resolutions of the singularity formed in the b = 0.59 flow. 
The size of the region where the surgery is performed is different in the two cases, and characterized 
by e = 0.07 (orange) and 0.14 (blue), but the resulting flows appear very similar after a short time. 
The geometries are plotted for intervals 0.002 in A after the resolution. Both flows eventually tend to 
flat space, although this is not shown here for clarity. 



of e therefore determines the size of the region where surgery is performed. This is smaller 
the closer one takes Aq to the singularity. In figure 14 we show two flows for b = 0.59 which 
use different values of Aq, giving e = 0.07 and 0.14. We note that both do indeed flow to 
flat space, and furthermore, the large scale form of the How appears to be insensitive to the 
details of the resolution. 
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